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ABSTRACT 

A five-parameter expression for a smoothly broken power law is presented. It is used to fit 
Gamma-Ray Burst spectra observed by BATSE. The function is compared to previously used four- 
parameter functions, such as a sharply broken power law and the Band et al. (1993) function. The 



presented function exists as a WINGSPAN routine at http://www. astro. su.se/^fclix/wing. htm 
KEYWORDS: Gamma-ray bursts. 

1. INTRODUCTION 

Studies of the spectral continuum of Gamma- Ray Bursts (GRBs) can give impor- 
tant clues for the understanding of the underlying emission processes responsible for 
the observed radiation. A number of studies have been done (Schaefer et al. 1992, 
1994, Band et al. 1993 , Ford et al. 1994). The burst emission is of a non-thermal 
character and most of the energy is released in the gamma-ray band (0.1-1 MeV). 
Spectra from missions prior to the Compton Gamma-Ray Observatory (CGRO) 
with its Burst and Transient Source Experiment (BATSE) were often fitted with an 
"optically thin thermal bremsstrahlung" photon spectrum, Ne(E) oc E~ 1 e~ E / E ° , 
combined with a power law at larger energies. Schaefer et al. also used a bro- 
ken power law in fitting BATSE bursts. Later, Band et al. showed that nearly 
all time-integrated BATSE spectra were well fitted by the four-parameter Band 
function consisting of two power laws, smoothly joined by an exponential roll-over. 
Below, we introduce a five parameter function, which consists of two power laws 
smoothly joined together. The sharpness of the turn-over is determined by the fifth 
parameter. 

2.1 SMOOTHLY BROKEN POWER LAW 

To model the time-integrated, background-subtracted photon count spectra, we 
want a function describing a smoothly broken power law, having a low energy 
power law with spectral index a, smoothly transferred into a high energy power law 
with spectral index f3. We follow the idea in Preece et al. (1996) and define the 
following properties: its logarithmic derivative varies smoothly from a to /3, with 
the transition described by a hyperbolic tangent function. 
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[Tabic 1.] The bursts sample. The bursts are identified by their date and trigger number (sec 
the Current BATSE Gamma- Ray Burst Catalogue, Meegan et al.). At and AEnergy are the 
chosen time and energy ranges used in the study. 



Date + 
Trigger No. 



LAD At 

(s) 



AEnergy 
(keV) 



Date + 
Trigger No. 



LAD At 

(s) 



AEnergy 
(keV) 



910430(# 130) 


6 


0-19 


29-1799 


910503(# 143) 


6 


0-10 


30-1766 


911127(# 1122) 


1 


0-38 


29-1774 


920210(# 1385) 


5 


0-49 


29-1988 


930506(# 2329) 


3 


0-14 


30-1756 



940217(# 2831) 





1-35 


33-1956 


950403(# 3491) 


3 


0-23 


30-1649 


970420(# 6198) 


4 


1-13 


27-1950 


980125(# 6581) 





45-59 


24-1980 


980203(# 6587) 


1 


0-34 


25-1855 



d\ogN E 
dlogE 



= £ tanh 



log(E/E ) 



(1) 



where Ne is the photon flux, E is the photon energy, Eq is the break energy, S is 
the width of the transition, £ = ((3 — a)/2, and (f> = ([3 + a)/2. Integrating (1) gives 



N E (E) = C 



E 



COSh ( ISiiM^ 



{<51nl0 



(2) 



where the width of the transition, in linear energy, is AE = E (10 s — 1). The 
normalisation is set by the constant C, which is the value of N E (E — E n ). The 
function is described by four parameters besides the normalisation; the two power 
law indices, a and /?, the break energy, E , and the energy range over which the 
spectrum changes from one power law to the other, 5 or AE. The advantage of this 
function over the Band function, for instance, which is described by four parameters, 
the normalisation, a, 0, and a parameter related to the break energy, is that the 
width of the transition region can be specified. 



3. GRB SPECTRA OBSERVED BY BATSE 

We have studied a sample of 10 strong GRBs, for which the Smoothly Broken 
Power Law (SBPL) above gives reasonable fits. They are used to illustrate the 
comparison between the Band model and the SBPL model. In Table 1, the set of 
bursts is displayed, together with the time and energy range used and the number 
of the Large Area Detector (LAD) used. The following three models were used: 
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[Tabic 2.] The results of the modelling of the bursts. A is the normalisation of the spectrum, 
Oi and (3 are the low and high energy power law photon indices, Eq is the break energy, and 
AE is the size of the turn-over spectral range. The goodness of fit is given by the % 2 -valuc 
together with the number of degrees of freedom, dof. See the text for details. 
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(1) A sharply Broken Power Law (BPL) consisting of two power laws connected 
at the break energy with three parameters besides the normalisation: the photon 
indices of the low and high energy power laws, a, j3, and the break energy, Eq. 

(2) The Band et al. (1993) function (Band) with three parameters besides the 
normalisation: a, Eq, (3. The spectrum below E — (a — (3)Eq is an exponentially 
cut off power law, N E (E) oc E a e~ E / Eo . A second power law, N-e(E) oc E@, is 
connected to this function above this energy. Note that -E pca k = (2 + a)E is 
usually used to characterise the break in the Band function. 

(3) The Smoothly Broken Power Law (SBPL) given by Equation (2) with its four 
parameters besides the normalisation. 

The results of the fitting are shown in Table 2. The general trends of the sample 
studied are the following. 

(1) The Band model gives the largest break energy (except for # 2831, which, 
however, has large errors). 

(2) The SBPL model has a somewhat larger break energy than the BPL model. 

(3) The low energy power law in the Band model is harder than the one in the 
SBPL model, which is marginally harder than the one in the BPL model. 

(4) The high energy power law in the BPL model is the hardest and in the Band 
model it is the softest. 

The main difference between the SBPL model and the Band model is the extra 
variable parameter describing the energy range of the changing of the power law 
index. As bursts can have an actual turn-over that is sharper than the one given by 
the fixed exponential turn-over in the Band function, which is determined solely by 
Eq, the resulting fits in Table 2 differ. The smaller the ratio between AE and Eq 
in the SBPL model is, the larger the difference is between Eq in the Band model 
and in the SBPL model. The best examples of this trend are # 6581 and # 143 
compared to # 1385 and # 1122. The smooth turn-over in the Band model allows, 
among other things, the low energy power law to be harder than what is found 
by the BPL model. However, extra parameters are not always advantageous. In 
some cases, not shown in Table 2, the AE parameter in the SBPL model cannot be 
constrained. 
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